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Introduction: The origin of geologically recent gul- 
lies on Mars has remained controversial since the discovery 
of these features by Malin and Edgett in 2000 [1]. Numer- 
ous models have been proposed which invoke various physi- 
cal processes as well as various agents of erosion to explain 
the origin of the Martian gullies. Hypotheses to explain the 
formation of the gullies invoke shallow liquid water aquifers 
[2], deep liquid water aquifers [3], melting ground ice [4, 5], 
snowmelt [6, 7, 8], dry landslides [9], and carbon dioxide 
aquifers [10]. We test the validity of such gully formation 
mechanisms by analyzing data from the Mars Global Sur- 
veyor and Mars Odyssey spacecrafts to uncover trends in the 
dimensional and physical properties of the gullies and their 
surrounding terrain. 

A similar study has previously been completed for gul- 
lies located in the southern hemisphere of Mars [11]. The 
work presented here focuses exclusively on gullies in the 
northern hemisphere based on the identification of 136 Mars 
Orbiter Camera (MOC) images containing clear evidence of 
gully landforms, distributed in the northern mid and high 
latitudes. These sites have been analyzed in combination 
with Mars Orbiter Laser Altimeter (MOLA), Thermal Emis- 
sion Spectrometer (TES), and Gamma Ray Spectrometer 
(GRS) data to provide quantitative measurements of numer- 
ous gully characteristics. Parameters measured include ap- 
parent source depth and distribution, vertical and horizontal 
dimensions, slopes, compass orientations, near-surface ice 
content, and factors controlling present-day climatic condi- 
tions. 

Methodology: Images from the Mars Orbiter Camera 
(MOC) have been systematically examined in search of gully 
features. Edgett et al. [12] report that the gullies nominally 
occur poleward of 30° in both the northern and southern 
hemispheres. Hence all images from the MOC narrow-angle 
camera between 30°N to 90°N from mission phases AB1 
through R15 (September 1997 through March 2004) were 
individually examined for the presence of gullies. Quantita- 
tive measurements of linear distances were then extracted 
from the MOC images, elevation data was extracted from the 
MOLA data, and TES measurements of thermal inertia were 
used in conjunction with modeled surface temperatures to 
derive subsurface temperature profiles. Gully locations were 
also compared with GRS data to investigate whether gully 
occurrence is correlated with the presence of near surface 
ground ice. 

Results: Based on this data, we present numerous 

trends with respect to the gully systems which must be ex- 
plained by any viable model of gully formation. Findings 


reported here generally support earlier results discussed by 
Heldmann and Mellon [11]. 

We find that the number of gully systems is greatest be- 
tween 30°N-45°N and tapers towards more poleward lati- 
tudes (Figure 1). Gullies form in a variety of terrain types as 
84% form in craters, 4% in knobby terrain, 3% in valleys, 
7% in other terrain types, and 2% in unknown terrain (no 
wide angle MOC context image). The geographic distribu- 
tion of gullies is shown in Figure 2. Northern hemisphere 
gullies are most concentrated in Utopia Planitia, Acidalia 
Planitia, Tempe Terra, and Arcadia Planitia. Additionally, 
gullies are found on all slope orientations at all latitudes but 
overall are preferentially found on equatorward facing slopes 
in the northern hemisphere (Figure 3). 

We also examine subsurface temperature and pressure 
fields to assess the possibility of liquid aquifers to feed the 
gullies. Using thermal conductivities derived from TES 
measurements as well as modeled surface temperatures, we 
calculate the depth to the 273K isotherm (liquid water melt- 
ing point) for all 136 locations of gully occurrence. We find 
the depth to the 273K isotherm is less than 200 meters for 
90% of the gully systems (Figure 4). Most gully alcove 
bases are below this required depth since based on measure- 
ments of 42 alcoves, the average depth to the alcove base is 
350 m. Using thermal conductivities derived from TES 
measurements as well as modeled surface temperatures, we 
find that ~93% of these gully alcove bases lie at depths 
where subsurface temperatures are greater than 273K and 
~7% of these alcove bases lie within the solid water regime. 
Interestingly, none of the gully alcoves lie within the tem- 
perature-pressure space of liquid C0 2 . Additionally, gullies 
generally are found in areas of low water content based on 
the GRS data (Figure 5). Although GRS is only sensitive to 
the upper meter of soil, this finding is consistent with the 
relatively low thermal conductivities as derived from TES 
data and the implication of a dry overburden layer which 
results in the relatively low (<200 m) depths to the 273K 
isotherm in these regions (Figure 4). 

No individual model yet explains all of the observed 
gully features, but based on this study and the similar analy- 
sis completed for the southern hemisphere [11] we can place 
additional constraints on future models and rule out several 
proposed hypotheses. Based on a comparison of such meas- 
ured gully features with predictions of the various models of 
gully formation, we find that the carbon dioxide, melting 
ground ice, dry landslide, deep aquifer and snowmelt models 
are the least likely mechanisms to describe the formation of 
the Martian gullies. Although some discrepancies still exist 
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between prediction and observation, the shallow aquifer 
model still remains as the most viable theory. 
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Figure 1. Histogram showing the number of MOC images 
containing clear evidence of gullies per 3° latitude bin in the 
northern hemisphere. 
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Figure 2. Geographic distribution of Martian gullies. 
Southern hemisphere gully locations are derived from 
Heldmann and Mellon (2004). 




Figure 3. The orientations of the Martian gully systems are 
plotted for three latitude bins as well as for 30°N-72°N. The 
labels N, S, E, W refer to the direction from which the gul- 
lies flow. Gullies plotted above the X axis flow from the 
north towards the equator. Gullies plotted below the X-axis 
flow from the south towards the pole. Contours are in inter- 
vals of two gullies within each angle bin of 1°. 



Figure 4. Calculated depth to 273K isotherm for northern 
hemisphere gully locations based on TES surface tempera- 
tures and thermal conductivities. 



Figure 5. Location of Martian gullies with respect to water 
content based on GRS data; percentages refer to water ice 
content. Southern hemisphere gully locations are derived 
from Heldmann and Mellon (2004). 




